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NOMENCLATURE 


a,b 


f 

F 

G 

h 


J 

k 


m 

m 


dimensions  of  computational  region,  Eqs.  (53),  (54) 
wave  velocity,  Eq.  (42) 

(oU*  -  uu)/(  a  -  w/c ) ,  Eq.  (25) 

w 

wave  propagation  speed  normalized  by  uT 

damping  coefficient,  Eqs.  (42),  (43),  (44),  and 
(47) 

unperturbed  channel  drag  coefficient 

pressure  drag  coefficient,  Eq.  (38) 

bending  rigidity,  Eq.  (46),  also  generalized 
vector,  Eq.  (58) 

spectral  density,  also  modulus  of  elasticity 
generalized  vector  in  Eq.  (58) 
a  function  defining  the  LES  filter,  Eq.  (71) 
velocity  vector,  Eq.  (58) 

filter  kernel  function,  Eqs.  (72)  and  (73) 

value  of  vertical  coordinate,  X2  at  the  flexible 
surface  (Fig.  1)  normalized  by  5  also  thickness 
of  wall  material 

Jacobian  of  coordinate  transformation 

spatial  wavenumber  in  Fourier  transform  space;  also 
spring  constant 

mass  per  unit  surface  area,  Eq.  (42) 
mass  flux 
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NOMENCLATURE  (continued) 


1 


*1»  M2 


Nx ,  Nz 


mass  terms,  Eqs.  (43)  and  (44) 

Mj/pT 

denotes  a  time  step  in  a  numerical  scheme 
membrane  tension  in  x  and  z  directions 
pressure  fluctuation  normalized  by  pu^ 
mean  pressure  normalized  by  pu^ 
resolvable  pressure  in  LES 
pressure  fluctuation  at  wall 
transformed  velocity,  Eq.  (84) 
u  tp/v 

two-point  correlation,  Eqs.  (50)  and  (51) 
velocity  gradients,  Eq.  (91) 


t  time  normalized  by  3/uT 

u^,i=l,2,3  velocity  component  (Fig.  1)  normalized  by  uT 


(u?)  ' 


<(u"2)>1/2 

1 

friction  velocity,  / tw/p  °f  flow  in  a  rigid  wall 
channel 

mass  average  velocity  of  channel  flow 


x^ , i  =  l ,  2 , 3 


Un/UT 

( 3<u, >/3x~ ) I 

1  1  'wall 

coordinate  (Fig.  1)  normalized  by  3 
the  vector  (x3,x2,x3) 


L.  .1  -  .-V  ^  -  f.  _  » 


NOMENCLATURE  (continued) 


coordinates  for  analysis  of  solid  wall  elastic 
behavior,  Eqs.  (42)  thru  (47)  and  associated  sketch 

distance  from  a  wall 

wave  number  in  x-^  direction  of  prescribed  wall 
shape  normalized  by  1/3 

wave  number  in  x-j  direction  of  prescribed  wall 
shape  normalized  by  1/3 

mean  channel  half-width 

Kronecker  delta 

ratio  of  change  in  flow  energy  to  energy  of  flow 
through  unperturbed  channel 

numerical  interval,  Eqs.  (74)  and  (75) 

antisymmetric  tensor 

displacement  of  wall  from  undisturbed  plane 
(Fig.  1)  normalized  by  3 

amplitude  of  prescribed  wall  displacement 
normalized  by  6 

direction  of  travel  of  wall  wave  relative  to  the 
xj  axis 

wavelength  of  wall  shape,  normalized  by  3 
fluid  viscosity 

1  -  X2 ,  Eqs.  (19)  thru  (22);  also  generalized 
coordinate,  Eqs.  (48)  and  (49),  Eqs.  (61),  etc. 

fluid  density 

[(w/c)2  -  i  (  wRe  )]1/2,  Eq.  (24)  • 
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NOMENCLATURE  (concluded) 


tw  fluid  shear  stress  at  wall 

t  transformed  time  variable,  Eq.  (61) 

a)  frequency  of  wall  motion  normalized  by  uT/$ 

Special  Notations 

<  >  denotes  an  average  over  an  horizontal  (xj ,x3)  plane 

for  LES  calculations  plus  an  average  over  time  for 
other  calculations 


(  ") 

R 

I 

3/  3x 
6/  6x 

A/  At 


denotes 

denotes 

denotes 

denotes 

denotes 

partial 

denotes 

partial 


a  periodically  varying  quantity 

a  real  part  of  a  complex  variable 

the  imaginary  part  of  a  complex  variable 

a  partial  spatial  derivative 

a  finite  difference  approximat ion  of  a 
derivative 

a  finite  difference  approximation  of  a 
time  derivative 


Superscripts 

'  denotes  the  nonperiodic  part  of  a  fluctuating 

turbulence  quantity 

"  denotes  a  fluctuating  turbulence  quantity 

+  denotes  a  length  in  "wall  units,"-  e.g.,  uTx/v 

*  denotes  a  dimensional  quantity 


Subscripts 

r,  i  denote  real  and  imaginary  parts  when  used  with  the 

quantities  A  and  a 
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INTRODUCTION 


Kramer  (references  1  and  2)  showed  that  the  application  of  a 
compliant  coating  to  a  wall  could  lead  to  a  considerable  reduction 
in  the  skin  friction.  He  focused  on  reducing  the  drag  by  delaying 
the  onset  of  transition  to  turbulent  flow.  Subsequent  attempts  to 
replicate  his  work  and  to  determine  the  causal  mechanisms  have  been 
only  partially  successful. 

A  number  of  studies  of  the  interaction  between  a  flow  and  a 
moving  wavy  surface  have  been  made;  these  have  been  principally 
motivated  by  the  air-sea  interaction  problem.  Benjamin  (reference 
3),  among  others,  showed  that  the  pressure  distribution  over  such  a 
surface  is  not  in  phase  with  the  shape  of  the  surface;  this  gives 
rise  to  a  momentum  transfer  between  the  fluid  and  the  surface,  the 
pressure  drag.  Although  these  studies  employed  an  inviscid  theory, 
they  have  played  an  important  role  in  the  studies  of  the  air-sea 
problem . 

Kendall  (reference  4)  investigated  the  turbulent  flow  over  a 
wavy  surface  experimentally  and  verified  some  of  the  predictions  of 
the  inviscid  theory  but  also  found  some  differences. 

Ff owcs-Wi 1 1 iams  (reference  5)  and  Blick  (reference  6)  studied 
the  possibility  of  drag  reduction  by  the  introduction  of  wavy  motion 
in  the  wall  by  means  of  a  theory  similar  to  the  one  that  will  be 
presented  in  this  report  and  came  to  the  conclusion  that  significant 
reductions  in  the  drag  should  indeed  be  possible. 


Hanratty  and  his  co-workers  at  the  University  of  Illinois  have 
conducted  a  number  of  studies  of  turbulent  flow  close  to  wavy  walls 
to  obtain  insight  into  the  influence  of  a  compliant  surface  on  a 
turbulent  flow  field.  Progress  has  been  made  in  calculations  of 
turbulent  flow  over  small-amplitude  stationary  waves  for  which  no 
flow  separation  occurs  (reference  7)  as  well  as  for  large-amplitude 
waves  (reference  8).  A  critical  issue  in  both  cases  is  the 
specification  of  the  wave  induced  variation  of  the  properties  of  the 
turbulence . 

The  work  reported  herein  was  motivated  by  the  possibility  of 
using  deliberately  introduced  wall  waviness  to  reduce  the  drag  on  a 
surface  underlying  a  turbulent  flow.  The  tool  used  in  this  study  is 
a  computational  one  —  large  eddy  simulation  (LES)  —  in  which  the 
large-scale  motions  of  the  flow  are  explicitly  computed  while  the 
mall-scale  motions  are  simulated  through  a  model.  The  method  has 
been  shown  to  capture  most  of  the  important  features  of  turbulent 
flow  over  flat  rigid  walls  including  the  structural  features  that 
are  responsible  for  most  of  the  turbulence  dynamics  in  the  near-wall 
region  (reference  9). 

In  this  report,  large  eddy  simulations  of  turbulent  flow  over  a 
prescribed  travelling  wavy  wall  will  be  presented.  The  simulations 
were  carried  out  with  linearized  boundary  conditions.  The  results 
indicate  that  significant  reduction  in  the  pressure  gradient  needed 
to  drive  the  flow  through  the  channel  are  indeed  possible.  However, 
additional  work  is  required  to  drive  the  wall  to  produce  this  effect 
and  the  effort  expended  is  a  considerable  fraction  of  the  reduction 
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in  the  work  needed  to  drive  the  through  flow.  A  linear  analysis 
explaining  the  principal  features  of  the  effects  observed  will  also 
be  presented. 

In  the  following  section,  the  principal  approximations  used  in 
this  work  will  be  introduced;  these  include  the  linearizations  of 
the  boundary  conditions  on  which  the  remainder  of  the  work  was 
based.  This  is  followed  by  two  sections  in  which  the  analysis  is 
presented.  The  first  gives  an  analysis  of  the  Stokes  layer  which 
was  found  to  exist  near  the  wall  by  Ffowcs-Williams  (reference  5) 
and  was  also  observed  in  the  LES  results.  This  is  followed  by  an 
analysis  of  the  effect  of  the  Stokes  layer  on  the  pressure  gradient 
and  an  estimate  of  the  work  done  by  the  wall.  Then  the  results  of 
the  large  eddy  simulations  are  presented  and  shown  to  be  in 
substantial  agreement  with  the  predictions  of  the  analysis. 

Finally,  conclusions  relating  to  the  possibility  of  drag  reduction 
are  discussed  and  directions  for  future  work  are  given. 

2.  BOUNDARY  CONDITIONS  FOR  FLEXIBLE  MOVING  WALLS 

2 . 1  Approximations 

Consider  flow  over  a  surface  which  is  able  to  stretch  and  move 
in  the  direction  normal  to  itself.  The  surface  is  described  by  the 
displacement  nfx^x^jt)  from  its  undisturbed  plane  shape,  where  x^ 
and  x-j  are  the  horizontal  coordinates  in  the  streamwise  and 
transverse  directions,  and  t  is  the  time.  The  exact  boundary 
conditions  for  the  viscous  flow  at  the  surface  are  that  the  fluid 


velocity  relative  to  the  surface  vanishes.  For  small  amplitude  of 
the  wall  motion,  a  convenient  approximation  is  to  expand  the 
velocity  boundary  condition  in  a  Taylor  series  about  the  undisturbed 
wall  plane  and  apply  the  boundary  conditions  on  that  plane.  The 
calculations  can  then  be  performed  in  a  channel  with  parallel  walls 
and  complicated  time-dependent  coordinate  transformations  can  be 
avoided . 

The  process  of  simplification  is  begun  by  approximating  the 
velocity  at  the  flexible  wall,  whose  undisturbed  position  is  *2  =  ^ ' 
by  two  terms  of  a  Taylor  series. 


3u  • 


ui(x1,h,x3,t)  =  ui(x1,l,x3,t)  +  n(xlfx3,t)  ) 


(1) 


x2  =  ! 


where  i=l,2,3  and  the  coordinates  and  velocity  components  are 
defined  in  Figure  1.  The  coordinates  are  normalized  by  the  mean 
channel  half-width. 

The  condition  of  zero  velocity  relative  to  the  surface  at 
X2  =  h  =  1  +  rif  yields  the  exact  boundary  condition: 


u ^ (x^  ,  h,x3  ,t )  =  0?  i  =  1  and  3 


(2) 


u2(x3 ,h,x3,t)  =  3n/3t 


(3) 


where  tangential  motion  due  to  stretching  has  been  neglected 


+1 


To  obtain  the  required  approximate  conditions  at  X2  = 
Equation  (1)  is  substituted  into  Equations  (2)  and  (3)  to  get: 


and 


3ut 

u  .  (x. ,l,x.,t)  =  -n  ( - ) 

3x2 


;  i  =  1  and  3 


(4) 


x2  =  1 


3u. 


u _ ( x 1 ,1 rx. ft)  =  —  -  n  ( — -) 

3 1  3x2 


(5) 


x2  =  l 


To  simplify  the  boundary  conditions,  the  velocity  is  considered  as  a 
combination  of  the  time  and  horizontal  mean  value  and  a  perturbation 
about  that  mean: 


u^ (Xj fx2,x^ ,t)  =  <u^(x2)>  +  u ^ (x j  ,x2  ,x^ f t )  ;  i  =  1,2,3 


(6) 


where  <  >  denotes  the  horizontal  and  time  average.  Then 


3u 

3x 


1 

2 


3<u  .  > 
1 


+ 


(7) 


It  is  assumed  that  n  <<  T  and 


0(n) 


(8) 


Therefore,  to  first  order  in  n 


/  3  <  U  i  >  v 

ux  ■  -  n(“3T7-). 

2  lx2=  1 


;  i  =  1,  3  (9) 


and 


u 


2 


3  n 

at 


(10) 


The  assumption  that  3u|/3x2  is  small  may  be  a  severe 
approximation  and  its  validity  should  be  checked  by  calculations 
which  do  not  use  it.  It  should  be  noted  that  experimental  data  for 
flow  over  flat  plates  show  that  the  rms  value  of 
3u[/3x2  (<( 3u|/3x2) 2>1/2)  is  about  25%  of  3<u1>/3x2 

(reference  10). 

The  non-zero  horizontal  velocity  component  at  the  wall  replaces 
the  actual  surface  by  one  which  is  being  stretched  and  compressed. 

It  is  important  to  determine  whether  such  a  condition  accurately 
models  the  wavy  wall.  As  the  flow  moves  along  the  actual  wall,  it 
is  displaced  vertically,  while  remaining  attached  to  the  wall.  The 
vertical  component  given  by  Equation  (10)  approximates  the  vertical 


motion  of  the  wall.  The  displacement  effect  is  approximated  by 
Equation  (9)  which  indicates  that  the  velocity  at  x2  =  1  is  the 
linear  extrapolation  (or  interpolation)  of  the  velocity  from 
x2  =  1  +  n  to  x2  =  1 .  However,  Equation  (9)  produces  locally 
reversed  flow  for  u^  where  (l-x2)  is  positive. 

2.2  PARAMETERS 

The  calculations  performed  in  this  work  used  the  boundary 
conditions  (9),  and  (10)  with  the  wall  motion  prescribed  to  be  a 
progressive  wave: 

n(x1,x3,t)  =  n0cos  (aXj+  8x3  -  wt)  (11) 

where  the  wavenumbers  a  and  8  and  the  frequency,  w  are  related  to 
wavelength,  X,  wave  speed,  c,  and  angle,  9,  by: 

a)  =  2ttc/X 
a  =  (2tt/X)  cos  9 
8  =  (  2 it/ X )  sin  0 


-7- 


Previous  studies,  e.g.,  Norris  and  Reynolds  (reference  11), 
Hanratty  et  al  (reference  12 ),  Thorsness  et  al  (reference  7),  and 
Kelleher  and  Balasubramanian  (reference  13)  have  shown  that  under 
certain  conditions  the  use  of  these  linearized  boundary  conditions 
may  be  a  poor  approximation.  However,  Norris  and  Reynolds'  work 
involved  wavelengths  and  amplitudes  much  larger  than  the  ones  to  be 
considered  here  and  the  other  work  involved  stationary  wavy  walls  so 
that  generalization  should  be  made  with  caution.  The  linearization 
of  the  boundary  conditions  is  the  most  significant  approximation  in 
this  work. 

The  results  of  reference  13  indicate  that,  for  calculations  of 
laminar  flow  over  a  stationary  wave,  using  linearized  boundary 
conditions  gives  a  larger  amplitude  in  the  surface  shear  stress 
varation  than  for  calculations  using  the  nonlinear  boundary 
conditions.  As  the  wave  amplitude  is  increased,  the  amplitude  of 
the  shear  stress  variation  also  increases.  The  linearized  boundary 
conditions  predict  regions  of  negative  surface  shear,  indicating 
separation,  for  amplitudes  for  which  the  nonlinear  boundary 
conditions  indicate  that  the  flow  is  still  attached.  At  the  same 
time,  the  use  of  the  linearized  boundary  conditions  appears  to  have 
little  effect  on  the  phase  angle  of  the  predicted  surface  shear 
variation.  These  results  are  only  qualitative  and  cautionary  with 
regard  to  the  present  study  since  the  effects  of  flow  turbulence  and 
wall  motion  are  not  accounted  for. 


The  question  of  the  validity  of  the  linearized  boundary 
conditions  deserves  further  attention.  In  deriving  equation  (9)  it 
was  assumed  that  3<u^>/3x2  is  of  zero  order  in  n  while 
3u^/3x2  is  of  first  order  so  that  n(3uj/3x2)  is  of  second 
order.  This  assumption  was  based  on  the  observation  that 
experimental  data  for  flat  plates  show  that  the  rms  value  of 
3uj/3x2  is  about  25%  of  3<u^>/3x2  .  The  relative  magnitude  of 
3u'j/3x2  cannot  be  determined  a  priori;  therefore,  the  assumption 
that  it  makes  a  second-order  contribution  to  the  boundary  condition 
can  be  evaluated  only  by  direct  comparison  with  a  solution  with 
exact,  nonlinear  boundary  conditions.  This  evaluation  is  essential 
but  could  not  be  carried  out  within  the  scope  of  the  present  work. 

3.  ANALYTICAL  STUDIES 

Preliminary  numerical  simulations  indicated  that  the  principal 
effects  of  the  wall  motion  are  confined  to  a  thin  region  near  the 
wall;  this  will  be  demonstrated  later.  On  this  basis,  it  was 
decided  to  develop  a  simplified  analysis  for  the  near-wall  region. 
This  analysis  was  a  valuable  aid  in  determining  values  of  the 
parameters  which  could  yield  effective  alteration  of  the  flow  and 
viscous  drag.  Although  the  study  began  with  numerical  simulations, 
it  is  useful  for  the  purposes  of  this  presentation  to  begin  with  the 


The  analysis  consists  of  two  parts.  First,  the  effect  of  the 
wall  motion  on  the  fluid  near  the  wall  was  found  to  produce  a 
Stokes-like  layer.  The  Stokes-like  solution  was  then  used  to 
estimate  the  change  in  the  mean  flow  and  the  drag  reduction.  An 
energy  balance  was  then  used  to  estimate  the  benefits  and  costs  of 
the  wavy  wall.  Finally,  the  results  of  the  analysis  were  used  to 
choose  cases  for  LES  calculations. 


3.1  The  Linearized  Stokes  Laver  Solution 
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with  the  wall  motion  given  by 


M 


n(x1,x3,t)  =  nQe 


i(ax1  +  6x3  -  ait) 


(15) 


The  actual  solution  is  the  real  part  of  the  result. 

For  the  purpose  of  this  analysis,  the  moving  wall  is  the  upper 
wall  of  the  channel  (X2=l)  and  the  lower  wall  is  considered  to  be 
at  Substitution  of  eqs.  (14)  and  (15)  into  eqs.  (12)  and  (13) 

and  the  boundary  conditions: 


u1(x1,l,x3,t)  =  -n  IT 


(16) 


u~ (x. , 1  ,x_  ,t  )  =  — 
1  1  J  3t 


(17) 


u3 (x1  ,1 ,x3  ,t )  =  0 


(18) 
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yields 


p(x  )  =  t)  cAe~ ^  c ^  ^ 

2  o 

G.  (x.  )  =  -n  (U*  +  -2£  A)e-°5+n  3£  Ae  "U/c)5 

1  2  O  w  O 

U)  U) 

u_(x?)  =  -in  (<*)  -  A)e-0?  -  in  Ae'(^c)? 

2  2  0  O 


u3(x2) 


il2l  A(e-(“/c,{  -  e-°5) 


where 


is  given  by 


E 


a 


The  velocity  u. 


i  ( ax  +  ex  -  tot) 

Ui (xi  ,x2 ,x3 ,t )  =  R[ui(x2)e  ] 


R[u.  (x2 )]  cos  (a\^+  Bx3-  tot) 


X  [ui <x2  ))  sin  <aX1+  ex3-  wt) 


(26) 


where  R  indicates  the  real  part  and  I  indicates  the  imaginary 
part.  The  Reynolds  stress  is  the  average  over  time  and  xj  and 
x^  of  the  product  and  will  be  shown  later  to  be  the  agent  of 

the  principal  effect  of  wall  motion  on  the  turbulent  flow.  From  the 
above  results 
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A  A 


<U 1^  2>  =  ~  (R(U1)R(U2)  +  I(u1)I(u2)] 


=  -  n2(ac  +  U')[A.e-2ar?  - 
«  o'-  w ;  1  1 


-(  o  +to/c  )  C 

(A. cosai ?+Arsinai 5)  e  ] 


(27) 
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and  A^, 


or  and  denote  real  and  imaginary  parts  of 


I 


where  Ar 

A  and  a,  respectively.  This  Reynolds  stress  is  plotted  in  Figure 
2.  Note  that  it  undergoes  a  damped  oscillation  and  in  the  region 
nearest  the  wall  is  opposite  in  sign  to  the  unperturbed  Reynolds 
stress . 


3 . 2  Numerical  Solutions  for  Nonlinear  Stokes  Layer 

In  order  to  evaluate  the  effect  of  neglecting  the  convective 
terms  of  the  Navier-Stokes  equations  in  the  Stokes  layer  analysis,  a 
numerical  solution  of  the  problem  was  obtained  using  a  computer  code 
for  unsteady  viscous  incompressible  flows  (reference  14).  The  code 
uses  a  factored  semi-implicit  solution  algorithm  that  is  explicit  in 
convective  terms,  and  implicit  in  viscous  terms.  The  method  is 
second-order  accurate  in  both  space  and  time.  For  the  present 
calculations,  the  code  solved  the  two-dimensional  incompressible 
Navier-Stokes  equations  in  a  channel  with  the  lower  wall  rigid,  and 
periodicity  assumed  for  the  inflow  and  outflow  boundaries.  The 
boundary  conditions,  Eqs.  (9)  -  (11)  were  applied  to  the  upper  wall 
of  the  channel;  calculations  were  performed  both  with  the  complete 
nonlinear  Navier-Stokes  equations  and  with  the  convection  terms 
omitted.  All  quantities  were  nondimensionalized  in  "wall  units," 
i.e.,  in  terms  of  friction  velocity,  uT  and  viscosity,  v.  The 
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A  =  Au  /v  =  200 
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T 

The  results  for  the  wall  motion-induced  Reynolds  stress  near 
the  moving  wall  are  shown  in  Figure  3  for  three  different 
computational  mesh  sizes  and  the  analytical  solution 
[equation  (27)].  The  symbols  are  located  at  the  computational  mesh 
points.  Symbols  +  and  x  represent  the  finest  mesh,  open  symbols 
the  coarsest  mesh.  Some  discrepancy  between  the  analytical  and 
numerical  solutions  results  from  truncation  error  in  the  finite 
difference  formulation.  However,  a  more  significant  effect  is 
produced  by  the  neglect  of  the  convective  terms  in  the  linear 
solutions.  The  nonlinearities  reduce  the  amplitude  of  the  response 
of  the  fluid  in  the  Stokes  layer.  It  is  important  to  note  that  in 
the  cases  shown  in  Figure  3,  the  results  are  entirely  due  to  the 
interaction  of  the  wall  motion  with  the  adjacent  viscous  fluid.  No 
turbulence  is  included  in  the  computation. 
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3 . 3  Estimate  of  Effect  of  Stokes  Layer  on  Mean  Pressure 
Gradient 


3.3.1  Mixing  Length  Analysis 

The  second  step  in  the  analysis  is  to  obtain  an  estimate  of  the 
drag  reduction.  The  velocity  and  pressure  in  the  channel  flow  with 
a  wavy  wall  can  be  decomposed  into  three  components  (reference  15). 


ui(xrt)  = 


<u.(x2)> 


U. (x,t) 


+  Uf (x,t)  ;  i  =  l  ,2,3 


l  — 


(  28) 


p(x,t)  =  P(x1)  +  p { x  ,  t )  +  p ' ( x  ,  t )  (29) 

where  <u^>  is  the  mean  velocity,  u^  is  the  periodic  part  of  the 
solution  induced  by  the  wall  motion,  and  u|  represents  the 
turbulent  fluctuations:  P,  p,  and  p'  have  similar  meanings.  It 
is  assumed  that  the  periodic  component  is  nothing  more  than  the 
Stokes  layer  solution  derived  in  the  previous  two  sections. 
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The  nondimensional  Navier-Stokes  equations,  averaged  over  time 
and  horizontal  planes  are  (reference  11): 


3x, 


(<u1u2  +  uju^>)  = 


dP 
dx . 


1  3‘ 


Re  .  2 
t  3x2 


<V 


(30) 


where  it  has  been  assumed  that  <uiu2  +  ui  u2>  negligible,  i»e,, 

that  the  turbulence  and  the  Stokes  Layer  are  uncorrelated.  Assuming 
the  <u jU 2>  is  not  affected  by  the  presence  of  <u1u2>,  the  mixing 
length  eddy  viscosity  model  with  Van  Driest  correction  (reference 
16),  can  be  used  to  yield. 


—  [(1  +  c  -2-  <u  >)  <u1>] 

3x2  3x2  3x2 


Re  (-J-  <u  u  >  +  ^-) 
T  3x2  dXj 


(31) 


where 


Re 


=  Min  {  (  K  )  2  ,  [  .  41  (  1- |  x_  |  )2]  [l-< 


-Re  ( 1- 1 x0 I ) / 26  2 

•  T  2  it 


(32) 


Where  K  is  a  constant  which  will  be  determined  subsequently. 
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Calculations  of  the  <u^u2>  profile  were  performed  with  the 
assumption  of  symmetry  at  the  channel  centerline  which  implies  that 
both  walls  are  wavy.  Integrating  eq.  (31)  in  x2  with  the  conditions 
that  e  and  <u^u2>  are  zero  at  x2  =  -1  gives 


> 


[1  -  4eR 


e  (<ulu2> 

T 


+  X 


3  P 
2  3x. 


)i1/2l 


(33) 


The  <ux>  profile  (Figure  4)  was  then  obtained  by  numerical 
integration  of  Eq.  (33)  using  a  fourth  order  Runge-Kutta  scheme. 
The  required  values  of  <u^u2>  were  obtained  from  Eq.  (27).  The 
mass  flow  rate 


ft  /_1  <u^^>dx2 


(34) 


was  then  obtained  by  trapezoidal  rule  integration  of  the  calculated 
<uj>  profile.  The  value  of  the  constant  K  in  Eq.  (32)  was 
determined  by  trial  and  error  so  that  the  mass  flow  rate  calculated 
from  Eqs.  (33)  and  (34)  with  Cu^u^  =  0  would  equal  the  mass  flow 
rate  used  in  the  LES  calculations  to  be  discussed  subsequently.  The 
value  determined  in  this  way  was 


K  =  0.085 
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3.3.2  Momentum  Balance 


The  effect  of  the  wall  motion  was  determined  from  the  momentum 
balance.  The  time  averaged  momentum  equation  was  integrated  over  a 
control  volume  bounded  by  the  channel  walls  and  periodic  in¬ 
flow/out-flow  boundaries.  The  first  step  of  the  process  is 
averaging  over  horizontal  planes,  which  yields  Eq.  (30). 

The  next  step  is  to  integrate  over  the  width  of  the  channel  to  get. 

1 

-1 

(35) 


<u!  u2  + 


-2  5e_  +  1 

dxl  Re 


3<u1  > 

T5T 


Since  both  the  periodic  and  the  turbulent  Reynolds  stress  are  zero 
at  the  channel  walls,  Eq.  (35)  yields,  for  symmetric  flow 


dp_  =  J_  9<U1>|  _  _i_ 

dx.  Re  3x„  ,  Re  w 


(36) 


It  is  noted  that  since  the  walls  are  horizontal  there  are  no 
horizontal  pressure  forces  in  the  momentum  balance.  Thus,  the  only 
effect  of  the  waviness  that  is  felt  by  the  flow  is  that  of  the  wall 
motion  on  the  mean  velocity  gradient  at  the  wall. 

In  order  to  determine  the  effect  of  wall  motion  on  the  velocity 


In  order  to  determine  the  effect  of  wall  motion  on  the  velocity 
gradient  at  the  wall,  Eq.  (33)  was  integrated  using  Eq.  (36)  and 
different  values  of  dp/dx^.  The  value  of  the  dp/dx^  was  adjusted  by 
Newton's  iteration  method  so  the  calculated  mass  flux  given  by  Eq. 
(34)  was  equal  to  the  mass  flux  of  the  unperturbed  channel.  The 
resulting  increment  in  dP/dx-^  is  assumed  to  be  twice  the  value  that 
would  be  obtained  with  one  moving  wall  and  one  rigid  wall. 

3.3.3  Results 

The  results  for  five  specific  cases  are  listed  in  Table  1.  The 
table  includes  results  from  numerical  integration  of  equations  (33), 
and  (34)  ,  and  from  LES  calculations  which  will  be  discussed 

subsequently.  The  values  listed  for  the  changes  in  dP/dxj,  from 
Eqs.  (33)  and  (34)  are  half  the  values  actually  calculated  in  order 
to  compare  with  the  LES  results  which  were  obtained  for  a  channel 
with  one  rigid  and  one  moving  wall. 

The  mixing  length  analysis  was  used  to  calculate  the  wall 
motion  effects  over  a  continuous  range  of  the  various  parameters.  A 
few  cases  were  then  selected  for  LES  calculations.  Case  1  is  the 
reference  case  for  which  the  parameters  were  listed  in  section 
3.2.  The  value  chosen  for  c  corresponds  approximately  to  the 
convection  velocity  of  large  scale  eddies  (reference  17).  The  value 
of  five  wall  units  for  n  +  is  roughly  the  thickness  of  the  viscous 
sublayer  and  was  considered  to  be  the  limit  of  applicability  of  the 
linearized  boundary  conditions.  The  wavelength  of  200  wall  units 


approximates  the  streamwise  spacing  of  large  eddies  near  the  wall. 

Case  2  is  the  same  as  Case  1  except  for  the  negative  value  of  c 
which  corresponds  to  a  wave  moving  upstream.  Case  3  corresponds  to 
wall  waves  rotated  90°  and  moving  in  the  spanwise  direction.  Case  4 
is  similar  to  case  2  to  determine  the  effects  of  the  magnitude  of 
the  wave  speed. 

The  effect  of  the  wall  motion  on  the  mean  velocity  profile  as 
calculated  by  the  integral  theory  is  shown  in  figure  4  for  case  2. 
The  figure  illustrates  that  the  effect  of  the  wall  motion  on  the 
mean  velocity  profile  is  felt  throughout  the  channel.  This  will  be 
discussed  further  when  the  LES  results  are  presented.  The  inset  in 
Figure  4  illustrates  that  the  3<U}>/3X2  at  the  wall  is  lower  for 
the  wavy  wall  than  for  the  rigid  wall.  The  momentum  balance  shows 
that  this  reduces  dP/dx^.  This  figure  also  indicates  the 
resolution  required  to  correctly  calculate  the  flow  by  LES;  the 
computational  mesh  should  have  a  point  0.2  wall  units  from  the  wall. 

An  important  observation  is  that  the  wall  motion  can 
significantly  affect  the  pressure  gradient.  The  simplified  theory 
presented  here  gives  a  larger  pressure  gradient  reduction  than 
LES.  A  major  reason  for  this  difference  is  that  the  Reynolds  stress 
is  overpredicted  by  the  linear  theory  (Fig.  3). 

The  decrease  in  the  driving  pressure  gradient  can  be  explained 
as  follows.  The  wall  motion  sets  up  a  Stokes-like  layer  in  its 
immediate  vicinity.  This  layer  has  a  "Reynolds  Stress" 

<u^U2>  opposite  in  sign  to  the  Reynolds  Stress  of  the  turbulence 
itself.  The  added  "Reynolds  Stress"  acts  to  decrease  the  turbulence 


production  near  the  wall  and  thereby  reduces  the  turbulent  Reynolds 
Stress  and  the  shear  force  on  the  wall.  Reduction  of  the  shear 
stress  must,  by  momentum  conservation,  be  accompanied  by  a  reduced 
pressure  gradient. 


3 . 4  Energy  Balance 

A  method  of  evaluating  the  net  effect  of  the  wall  motion  is  to 
examine  the  energy  balance.  The  mean  energy  equation  is  derived  by 
taking  the  scalar  product  of  the  momentum  equation  with  the  velocity 
and  averaging  over  time  and  a  control  volume  bounded  by  the  channel 
walls  and  periodic  inflow,  outflow  and  side  planes.  The  resulting 
energy  equation  is 
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(37) 


where  V  is  the  control  volume.  In  Eq.  (37)  <  >  indicates  time  and 
horizontal  plane  averages  evaluated  at  the  moving  wall. 

In  the  rigid-walled  channel,  the  energy  required  to  move  the 
fluid  through  the  channel  is  balanced  by  the  dissipation.  The 
moving  wall  contributes  to  the  energy  balance  through  the  work  done 
by  viscous  forces  and  by  the  fluctuating  pressure  at  the  wall  and  by 
changing  the  dissipation.  Thus,  any  reduction  in  the  required 
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energy  to  drive  the  flow  may  be  offset  by  the  work  required  to 
generate  the  Stokes  layer. 

Negative  values  of  the  wall  pressure-work  and  positive  values 
of  the  viscous-work  terms  represent  energy  added  to  the  fluid  by  the 
wall.  These  must  be  considered  penalties  to  be  deducted  from  the 
reduced  through-flow  pressure-work  requirement. 

The  various  components  of  the  energy  balance  are  shown  as 
functions  of  the  wave  speed  in  Figure  5.  The  energy  throughput 
increment  is  the  mean  velocity,  Um  times  the  change  in  the  mean 
pressure  gradient,  dP/dxj  from  its  unperturbed  channel  (half  the 
change  in  dP/dx-^  found  from  the  mixing-length  analysis  described  in 
section  3.3).  The  an<3  viscous  wall  work  terms  are  evaluated 

using  the  linear  theory.  The  change  in  dissipation  was  not 
computed;  it  is  estimated  that  the  dissipation  is  reduced  somewhat 
by  the  wavy  wall  so  the  results  should  be  conservative.  Clearly,  a 
large  part  of  the  change  in  the  work  required  to  drive  the  through 
flow  is  offset  by  the  viscous  work  done  by  the  wall  motion.  This  is 
independent  of  the  direction  of  the  motion.  On  the  other  hand,  the 
pressure  work  is  direction  dependent.  For  downstream  wave  motion  at 
moderate  wave  speeds,  the  pressure  work  is  negative,  indicating  that 
the  wall  absorbs  energy  from  the  fluid.  Thus,  the  fluid  drives  the 
wall  motion.  For  upstream  wave  motion,  the  wall  adds  energy  to  the 
fluid  through  the  <pu2>  interaction.  This  suggests  that  upstream 
wave  motion  is  not  possible  for  compliant  walls  because  it  cannot  be 
driven  by  the  flow. 
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The  net  effect  of  the  Stokes  layer  on  the  energy  required  to 
drive  the  flow  through  the  channel  is  shown  in  Figure  6.  A  small 
net  decrease  in  energy  required  is  produced  by  a  wall  moving  either 
upstream  or  downstream  at  moderate  speeds;  note  that  the  change  in 
dissipation  is  not  included.  The  theoretical  curve  shown  in  Figure 
6  displays  singular  behavior  around  c.  =  0.  This  is  an  inadequacy 
of  the  linearized  theory. 

The  results  shown  in  Figures  5  and  6  differ  from  results  given 
in  an  earlier  interim  report  (Ref.  18).  The  reason  is  that  in  the 
earlier  report,  the  <UjU2>  contribution  to  Eq.  (33)  was  reduced  by 
an  arbitrary  factor  of  0.5  to  account  approximately  for 
nonlinearities  (see  Figure  3).  The  arbitrary  reduction  was  not  done 
here.  Also,  in  the  previous  work  the  value  of  used  in  all 
calculations  of  quantities  from  the  linear  theory  was  the  rigid-wall 
value.  Since  the  slope  of  the  mean  velocity  profile  at  the  wall 
changes  along  with  the  mean  pressure  gradient,  the  present 
calculations  also  include  the  correct  value  of  U^,  given  by  Eq. 

(36).  Although  the  approximations  used  lead  to  results  that  are 
only  qualitative,  the  results  suggest  that  (a)  the  effect  of  a  wavy 
wall  is  to  reduce  the  energy  required  to  drive  the  flow  and  (b)  the 
reduction  is  small. 

In  addition  to  the  approximations  already  discussed,  the 
turbulent  pressure  fluctuation  was  not  included  in  calculating  the 
pressure  work  term.  The  effect  of  including  these  fluctuations  is 


unknown  as  the  theory  does  not  predict  the  change  in  the  turbulence 
quantities. 

As  the  linear  theory  predicts  a  larger  magnitude  of  reversed 

Reynolds  stress  near  the  wavy  wall  than  the  nonlinear  theory,  it  is 

expected  that  the  linear  theory  would  overpredict  the  reduction  in 

the  mean  pressure  gradient  and  the  viscous  work  and  the  pressure 

work  at  the  wavy  wall.  A  further  illustration  of  the  effects  of  the 

nonlinearities  can  be  seen  in  the  variation  of  the  intensity  of  the 

velocity  fluctuations  produced  in  the  Stokes  layer.  Figure  7  shows 

~2 

the  variation  of  the  rms  velocity  fluctuations  /<u^>  near  the  wall 

for  the  linear  and  the  nonlinear  Stokes  layer  solutions  for  case 

~2 

2.  The  linear  solution  has  larger  (negative)  slope  (3/<  u^>/3x2) 
at  the  wall  than  does  the  nonlinear  solution.  This  suggests  that 
the  actual  viscous  work  would  be  smaller  than  predicted  by  the 
linear  theory.  However,  the  errors  in  the  predicted  quantities  are 
expected  to  be  small  and  to  follow  the  correct  trends.  Therefore, 
the  linear  theory  is  believed  to  be  a  useful  guide  for  LES 
calculations. 

This  study  indicates  that,  in  general,  the  wall  must  supply 
energy  to  produce  the  conditions  that  result  in  reduced  mean 
pressure  gradient.  The  question  of  whether  a  significant  net 
reduction  in  the  energy  required  to  drive  the  flow  is  possible  is 
the  essential  one.  This  investigation  indicates  that  the  amount  of 
energy  required  from  the  wall  is  roughly  equal  to  the  amount  of 
energy  saved  by  the  drag  reduction. 


To  make  this  calculation  possible  the  following  elements  are 
needed : 

1.  The  spectrum  in  space  and  time  of  the  pressure 
fluctuations  on  a  wall  in  natural  turbulent  flow.  This 
could  be  obtained  from  either  an  experiment  or  a  full  or 
large  eddy  simulation.  To  date,  only  limited  portions  of 
the  necessary  data  have  been  published,  namely  the  time 
spectra  of  the  pressure  fluctuations. 

2.  A  dynamic  linearized  model  of  the  wall  in  which  the  input 
is  the  pressure  fluctuations  on  the  surface  and  the  output 
is  the  surface  displacement.  This  is  just  the  dispersion 
relation  for  the  surface  and  needs  to  be  specified  from 
the  dynamics  of  deformable  bodies.  A  more  sophisticated 
model  might  include  the  effect  of  the  surface  shear 
stresses;  however,  these  are  not  likely  to  be  important. 

3.  A  model  for  the  effects  of  the  wall  surface  displacement 
on  the  flow.  The  Stokes  layer  model  presented  above 
should  be  adequate  for  this  task. 

The  suggested  combination  of  these  elements  would  produce  a 
model  capable  of  predicting  the  drag  in  a  compliant  wall  channel 
containing  a  turbulent  flow.  On  the  basis  of  the  experience 


described  in  this  report,  it  is  anticipated  that  this  theory  will  be 
useful  in  estimating  the  change  in  drag  and  the  properties  of  the 
surface  which  would  produce  the  greatest  drag  but  will  not  be 
quantitatively  accurate.  Unfortunately,  time  restrictions 
prohibited  the  present  study  from  investigating  this  model. 

3 . 6  Pressure  Drag 

The  momentum  balance  shows  that,  under  the  approximations 
employed,  there  is  no  pressure  drag  on  the  wall.  However,  an 
estimate  of  the  pressure  drag  on  the  actual  wavy  wall  can  be 
obtained  by  assuming  that  the  Stokes  layer  pressure  at  the  displaced 
wall  would  obtain  at  the  wavy  wall.  A  "drag  coefficient"  can  then 
be  found  by  integration. 
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Using  equations  (11)  and  (19),  this  becomes 
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The  viscous  drag  of  the  unperturbed  channel  is  in  equilibrium  with 
the  pressure  drop: 
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(40) 


so  that 


CD  /CD  =  -”o  aCV4  (41) 

P 

For  cases  1,  2,  and  4  of  Table  1,  the  pressure  drag  from  Eq.  (41)  is 
3.9  percent,  10.4  percent,  and  13.2  percent  respectively.  For 
spanwise  waves,  no  pressure  drag  is  produced.  The  LES  results  shown 
in  Table  1  will  be  discussed  subsequently. 

4.  COMPLIANT  WALL  DYNAMICS 

The  work  on  the  effects  of  moving  walls  with  prescribed  shapes 
on  turbulent  flow  provided  considerable  understanding  of  the 
interaction  between  wall  motion  and  the  fluid  flow.  The  principal 
purpose  of  this  work  was  directed  at  understanding  the  interaction 
of  a  compliant  wall  with  the  flow.  In  this  case,  the  equations  for 
the  wall  motion  must  be  solved  simultaneously  with  the  fluid 
equations. 
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If  only  streamwise  variations  in  wall  shape  are  allowed,  the 
dynamic  equation  satisfied  by  the  displacement  y  =  n(x,t)  of  the 
plane  surface  bounding  an  elastic  medium,  including  linear  damping, 
is  the  second  order  equation  (Benjamin,  reference  3) 
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velocity  of  free  waves 
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Equation  (42)  can  easily  be  extended  to  allow  variation  in  the 
spanwise  direction. 

Yang  and  Heller  (reference  19)  analyzed  a  viscoelastic  plate 
sandwiched  between  a  rigid  plate  and  flexible  membrane.  They 
obtained  a  pair  of  equations  similar  to  equation  (42)  containing 
several  additional  terms;  they  allow  for  variation  of  the  surface 
displacement  in  two  directions.  The  major  omissions  in  both  of 
these  approaches  are  the  shearing  stresses  and  rotary  inertia  in  the 
dispersive  material. 
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Accurate  accounting  for  all  possible  effects  is  shown  by 
Achenbach  (reference  20)  to  produce  a  set  of  three  rather 
complicated  equations  which  admit  many  modes  of  vibration,  each 
having  its  own  dispersion  curve.  However,  it  is  not  known  whether 
all  modes  are  important  in  the  fluid  interaction  problem.  The  modes 
associated  with  shearing  and  rotational  motions  are  expected  to 
interact  much  more  weakly  with  the  flow  than  the  modes  due  to 
bending  and  thickness  deformations  produced  by  a  fluctuating 
pressure . 

Under  the  compliant  Coating  Drag  Reduction  Program  sponsored  by 
the  U.S.  Office  of  Naval  Research,  a  considerable  amount  of  work  has 
been  done  to  determine  the  characteristics  of  elastic  surfaces  which 
could  interact  strongly  with  the  fluid  flow.  Much  of  the  effort  has 
been  directed  toward  development  of  mathematical  models  of  various 
combinations  of  layers  of  viscoelastic  materials.  Duncan  and  Hsu 
(reference  21)  studied  single-and  double-layer  viscoelastic 
coatings,  and  estimated  the  material  properties  required  for  maximum 
interaction  with  pressure  fluctuations  of  the  flow.  Kalnins  and 
Evrensel  (reference  22)  demonstrated  how  a  composite,  compliant 
coating  could  be  modeled  by  a  one  degree-of -freedom  system,  similar 
to  a  membrane.  At  the  time  of  the  present  work,  a  definitive 
description  of  a  compliant  surface  had  not  been  developed.  The 
viscoelastic  sandwich  structure  studied  by  Yang  and  Heller  provided 
a  simple  means  of  testing  the  feasibility  of  coupling  the  LES  with  a 
compliant  surface.  The  equations  derived  for  the  surface  deflection 


due  to  an  imposed  pressure  distribution  represent  in  an  approximate 
way  surfaces  ranging  from  simple  membranes  to  viscoelastic  layers. 


The  compliant  wall  is  assumed  to  be  a  laminated  structure 
composed  of  a  layer  of  a  viscoelastic  material  between  a  stiff  plate 
and  a  thin  flexible  plate. 


The  base  plate  is  assumed  to  be  subjected  to  small  deflection,  while 
the  upper  plate  may  undergo  large  deflections  and,  in  addition,  may 
be  subjected  to  in-plane  forces.  According  to  Yang  and  Heller, 
experiments  have  shown  that  the  Poisson's  ratio,  v,  of  PVC  foam 
under  transverse  loads  is  nearly  zero.  This  allows  the  assumption 
that  the  middle  layer  can  be  treated  as  a  viscoelastic  spring.  The 
equations  of  motion  of  the  composite  plr-'  then  become 
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where  =  -j  p^h^  +  =  -j  p^h^  +  p2^2'  an<^  pl'  p2'  ^1'  ^2' 

n^(x,z,t)  and  n2(x,z,t)  are  the  masses  per  unit  volume,  thicknesses 
and  transverse  displacements  of  the  upper  plates  and  the  base  plate 
respectively,  and  p^  and  h-j  are  the  mass  per  unit  volume  and 
thickness  of  the  viscoelastic  medium.  The  quantity  C3  is  the 
damping  coefficient  of  the  viscoelastic  middle  layer. 


The  spring  constant  1<3  is  related  to  the  elastic  modulus  and  the 
thickness  of  the  middle  layer  by 
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(45) 


The  quantitities  and  D2  are  the  bending  rigidities  of  the 

upper  plate  and  the  base  plate,  respectively,  with  elastic  moduli 
and  E2  and  Poisson's  ratios  and  V2: 

Eihi 

D.  =  - ;  i  -  1,2  (46) 

1  1 2  (  1-  vf  ) 

The  function  p(x,z,t)  in  equation  (43)  is  the  boundary  layer 
pressure  fluctuation. 

If  the  base  plate  is  completely  rigid  so  that  D2  +  ®,  the  base 
plate  deflection  can  be  neglected,  eliminating  equation  (44). 
Equation  (43)  then  reduces  to 


Mi  L2 +  Cs +  k3'1  +  Di''4’’ 


2  2 

-Nx  “ — ?  ”  N  — — 5  =  p(x,z,t) 

X  8  xZ  Z  8zZ 


(47) 


In  what  follows,  all  lengths  are  nondimensionalized  with  the 
mean  channel  half-width,  "6  ,  densities  with  the  fluid  density,  p. 


and 


velocities  with  the  fluid  friction  velocity,  u  -  /  t  /p , 

T  W 

pressure  and  stresses  with  p  \x\.  The  values  used  for  the 
dimensional  quantities  are  listed  in  table  2. 

With  the  compliant  wall  displacement,  n  ,  included  as  a 
dependent  variable  along  with  the  fluid  velocity  components  and  the 
pressure,  equation  (47)  is  solved  simultaneously  with  the  channel 
flow  equations  by  writing  it  as  two  first  order  equations 
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which  are  solved  simultaneously  with  the  fluid  flow  equations  by  the 
(implicit)  Crank-Nicolson  method. 


5.  LARGE-EDDY  SIMULATION 

The  large-eddy  simulation  (LES)  code  used  in  this  work  is 
described  in  Reference  9.  The  equations  for  the  large-scale  flow 
field  are  obtained  by  integration  of  the  filtered,  three- 
dimensional,  time-dependent  Navier-Stokes  equations.  The  small- 
scale  field  motions  are  simulated  through  an  eddy-viscosity  model. 
The  program  permits  simulation  of  a  turbulent  channel  flow  at 
moderate  Reynolds  number  (Re^  = 


14000  )  . 


The  LES  code  was  shown  to  capture  most  of  the  important 
features  of  the  turbulent  channel  flow.  Calculated  results  are  in 
good  agreement  with  experimental  data  for  statistical  quantities. 
Low-  and  high-speed  streaks  alternating  in  the  spanwise  direction 
found  in  the  numerical  results  were  in  good  qualitative  agreement 
with  experimental  observations.  The  calculated  production  of 
turbulence  was  found  to  be  intermittent  in  a  manner  that  strongly 
resembles  that  seen  in  the  laboratory. 

The  impetus  for  the  present  work  was  the  capability  of  the 
three-dimensional  time-dependent,  numerical  simulation  to  provide 
detailed,  instantaneous  information  about  the  flow  at  many  spatial 
locations.  This  information  can  be  used  effectively  to  study  the 
structure  and  statistical  properties  of  the  flow.  For  the  case  of 
time-  and  spatially-varying  boundary  conditions,  this  capability 
provides  a  unique  opportunity  to  study  the  mechanics  of  the 
interaction  between  the  wall  motion  and  the  turbulent  flow. 

The  code  originally  was  developed  for  use  on  the  ILLIAC  IV 
computer  at  NASA/Ames  Research  Center.  Some  of  the  calculations  for 
prescribed  wall  motion  discussed  in  this  report  were  performed  on 
tha  .  computer.  The  code  was  modified  for  use  on  the  NASA/Ames  CRAY 
computer  where  the  remainder  of  the  calculations  were  done. 


The  LES  calculations  were  made  in  three  stages.  First,  a 
preliminary  calculation  was  made  to  evaluate  whether  the  choice  of 
wall  parameters  would  produce  observable  effects.  Then,  after  the 
analytical  studies  provided  guidelines  for  the  choice  of  wall 
parameters,  a  second  set  of  calculations  was  performed.  Both  stages 
employed  a  coordinate  mesh  that  consisted  of  64  points  in  the  x^  and 
X2  directions  (streamwise  and  normal  to  the  walls)  and  128  points  in 
the  X3  (spanwise)  direction.  In  stage  two,  both  prescribed  wall 
motion  and  compliant  wall  motion  were  studied.  In  the  final  stage, 
102  points  were  used  in  the  X2  direction  to  better  resolve  the  near 
wall  region.  Also,  the  terms  of  the  global  energy  balance  equation 
and  pressure  drag  were  calculated. 

In  the  following  sections,  the  results  of  both  the  coarse  mesh 
and  the  fine  mesh  LES  calculations  will  be  used  to  illustrate  the 
effects  of  the  wall  motion.  Some  of  the  results  of  the  coarse-mesh 
calculations  were  presented  in  reference  18. 


5.1  LES  Results  for  Prescribed  Wall  Motion 


The  four  cases  discussed  previously  (section  3  and  Table  1) 
were  computed  using  the  LES  code.  Since  the  sign  of  Reynolds 


stress  -<U3U2>,  where  u'J  =  uj_  -  <u^>,  for  a  steady 


turbulent 


channel  flow  is  negative  near  the  upper  wall,  wall  parameters  were 
selected  for  which  the  linear  theory,  Eq.  (27)  gave  positive  values 
of  -<u^U2>  ;  this  should  reduce  the  production  of  turbulence  and 
yield  drag  reduction. 
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The  first  computational  mesh  used  in  the  LES  code  was 
approximately  the  same  near  the  wall  as  the  coarsest  mesh  shown  in 
Figure  3.  The  finest  mesh  used  for  LES  calculations  had  102  points 
in  the  x2  direction,  and  near  the  wavy  wall  approximately 
corresponded  to  the  finest  mesh  shown  in  Figure  3.  The  details  of 
the  computational  mesh  are  listed  in  Table  3.  Note  that  the 
distance  from  each  wall  in  wall  units  is  included  in  the  table. 

These  values  are  only  approximate  since  they  are  based  on  the 
friction  velocity  of  a  rigid  channel.  For  the  channel  with 
different  boundary  conditions  on  opposite  walls,  the  mesh  spacing  in 
wall  units  will  differ  slightly  for  each  wall. 

Starting  from  an  initial  velocity  field  representing  the  flow 
in  the  rigid-walled  channel,  the  governing  equations  were  integrated 
forward  in  time  with  the  prescribed  moving  wall  conditions  until  the 
numerical  solutions  reached  statistically  steady  states.  The  mass 
flux  was  held  fixed  while  the  mean  pressure  gradient  was  allowed  to 
float . 

5.1.1  Effect  of  Wall  Motion  on  Pressure  Gradient 

In  Figure  8a,  the  pressure  gradients  for  the  rigid-walled 
channel  and  for  the  wavy-walled  channels  of  Cases  1,  2,  and  3,  of 
Table  1  are  shown;  the  coarse  mesh  was  used  for  these  cases.  Using 
the  fine  mesh,  cases  1  and  2  were  recalculated  and  an  additional 
case  (case  4)  with  a  higher-speed  upstream-traveling  wall  wave  was 
run;  these  results  are  given  in  Figure  8b.  The  terms  of  the  energy 
balance  equation  were  also  computed  in  the  fine  mesh  calculations. 


All  cases  have  an  initial  transient  period  wherein  the  solution 
adjusts  from  the  rigid  wall  solution  to  the  moving  wall  solution. 
Case  4  was  started  using  the  solution  from  the  last  time  step  of 
Case  2.  This  resulted  in  smaller  initial  oscillation  and  faster 
damping . 

The  coarse  and  fine  mesh  results  are  similar;  the  pressure 
gradient  reductions  agree  quite  well.  The  coarse  mesh  produces  a 
reduction  in  dp/dx^  of  about  5  percent  for  Case  1  and  nearly  20 
percent  for  Case  2,  while  the  fine  mesh  produces  reductions  of  about 
3.5  and  17  percent,  respectively.  For  case  3,  the  coarse  mesh 
calculation  produces  a  10  percent  reduction  in  dp/dxj.  The  high 
speed  case  4  calculated  with  the  fine  mesh  yields  a  greater 
reduction  (20  percent).  The  trends  are  consistent  with  the  linear 
theory  predictions. 

5.1.2  Effect  of  Wall  Motion  on  Turbulence 

Examples  of  the  turbulence  statistics  from  LES  calculations  are 
shown  in  Figures  9-16.  All  the  figures  show  results  of  the  fine  mesh 
calculations.  All  cases  show  similar  effects  of  the  moving  wall  on 
the  turbulence  statistics. 

In  Figure  9,  the  mean  velocity  profile  is  shown  for  Case  2 
using  the  fine  mesh;  the  behavior  is  similar  to  that  of  the  Stokes 
layer  theory  shown  in  Figure  4.  Of  particular  interest  is  the 
expanded  inset  which  shows  the  high  resolution  of  the  fine  mesh 
which  accurately  predicts  the  lower  value  of  the  d  <  Uj  >/dx2  at  the 
moving  wall. 


Figure  8  shows  only  400  time  steps  for  the  coarse  mesh  cases, 
1200  steps  for  the  fine  mesh  cases  1  and  2  and  900  steps  for  case 
4. The  total  run  length  for  the  coarse  mesh  cases  was  1600  steps;  the 
profiles  of  the  various  quantities  shown  in  Figures  9-16  were 
calculated  by  averaging  over  horizontal  planes  and  the  last  200  time 
steps.  Profiles  from  the  fine  mesh  calculations  represent  averages 
over  the  last  300  steps  after  a  total  run  length  of  1200  steps  for 
cases  1  and  2  and  900  steps  for  case  4.  Because  of  the  short  time¬ 
averaging  period,  the  profiles  suffer  from  some  statistical  noise. 
Nevertheless,  certain  trends  can  be  discerned  which  indicate  the 
effects  of  the  wall  motion  on  the  turbulence. 

In  Figure  10,  the  horizontally  averaged  turbulence  intensity 
11 2  1/2 

components,  Cu.^  >  '  from  Case  1  for  the  lower  and  upper  halves  of 
the  wavy-walled  channel,  respectively,  are  shown.  The  only 
significant  difference  between  the  two  halves  of  the  flow  occurs 
very  near  the  upper  (wavy)  wall,  where  the  u-^  and  U2  intensities 
increase  significantly.  The  effect  is  shown  more  clearly  in  Figure 
11,  where  the  turbulence  intensity  <uJ^>l/2  is  compared  for  the 
rigid  and  moving  walls  for  cases  1,  2,  and  4.  The  intensity  is  very 
high  at  the  moving  wall  in  all  cases  but  rapidly  falls  to  a 
variation  similar  to  that  near  the  rigid  wall. 

In  Figure  12  the  turbulent  shear  stress  for  Cases  1,  2,  and  4 
is  compared  for  the  rigid  and  moving  walls.  Near  the  wavy  wall  the 
shear  stress  changes  sign  as  expected.  Near  the  channel  center, 
statistical  noise  obscures  any  effects  that  might  be  present. 


In  Section  3  of  this  report,  it  was  shown  (Fig.  3)  that  the 
nonlinear  terms  of  the  Navier-Stokes  equations  have  a  large  role  in 
determining  the  effect  of  wall  motion  on  the  flow  Reynolds  stress 
and,  consequently,  on  the  drag  of  the  channel.  In  Figure  13,  the 
Reynolds  stress  from  the  nonlinear  Stokes  layer  solutions  for  cases 
1,  2,  and  4  is  compared  with  the  LES  results.  Clearly,  turbulence 
has  an  appreciable  effect  on  the  moving-wall/flow  interaction.  In 
particular,  for  the  downstream-moving  wall,  LES  produces  much  less 
<u^u^>  reversal.  Also,  it  is  clear  why  the  fine  mesh  LES  solution 
yields  less  reduction  in  the  channel  pressure  gradient  than  the 
coarse  mesh  solution.  The  high  resolution  mesh,  gives  a 
significantly  smaller  Stokes  layer  contribution  to  Cu'^u’^  . 

Two  point  correlation  functions 
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where  f"  represents  the  fluctuating  velocities,  u^"  for  i  =  1,2,3, 
or  the  pressure  fluctuation,  p  for  i  =  4,  are  plotted  in  Figure  14 
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at  two  vertical  locations,  one  near  the  lower  rigid  wall  and  one 
near  the  upper  wavy  wall  for  case  l.The  channel  centerline  is 
located  640  wall  units  from  each  wall.  In  this  section,  the  <  > 
indicates  an  horizontal  average.  Since  the  statistics  were  not 
time-averaged,  they  suffer  from  statistical  fluctuations, 
particularly  at  large  separation  distance.  The  longitudinal 
correlation  near  the  rigid  wall  shown  in  Figure  14(a)  is  like  that 
of  Reference  9  for  the  rigid-walled  channel.  For  small  separation 
distances,  the  correlation  for  u’j*  is  larger  than  that  for  the 
corresponding  transverse  components,  u^  and  U3  and  the  pressure,  and 
it  extends  over  much  longer  distances  than  do  the  other 
correlations.  This  is  also  true  in  the  upper  half  of  the  channel 
except  very  near  the  boundary.  As  the  wavy  wall  is  approached 
[Figure  14(b)],  the  correlation  becomes  dominated  by  the  wall 
motion.  At  x^=  .914,  ( y+  =  55  for  a  rigid  walled  channel)  some 

effect  of  the  boundary  motion  is  apparent  in  u'^  and  p.  A  smaller 
effect  is  discernible  in  the  u!|  ,  while  u!j  is  completely 
unaffected.  These  comparisons  show  that  the  effects  of  the  wall 
motion  are  felt  for  a  considerable  distance  into  the  flow. 

In  the  spanwise  direction,  the  two-point  correlations  show  very 
little  wall  motion  effect  except  very  near  the  upper  boundary  where 
the  uj  ,u£  and  p  become  highly  correlated  across  the  entire  span 
of  the  channel  due  to  the  two-dimens ional  nature  of  the  wall  motion. 

Energy  spectra  of  the  turbulence  are  shown  in  Figure  15  for  the 
same  two  locations  as  for  the  two-point  correlations.  The  effects 


of  the  wall  motion  are  evident  in  the  spectra  of  the  longitudinal 
and  vertical  energy  components  [Figure  15(b)]  where  sharp  peaks 
occur  at  the  wavenumber  of  the  prescribed  wall  motion. The  wall 
motion  makes  a  negligible  contribution  to  the  spanwise  energy 
spectra . 

The  same  comments  generally  apply  to  the  results  of  case  3 
except  that  the  largest  effects  occur  in  spanwise  two-point 
correlations  and  energy  spectra. 

5.1.3  Energy  Balance 

As  shown  in  Section  3.4,  the  presence  of  a  moving  wavy  wall 
introduces  two  new  terms  to  the  kinetic  energy  integral  equation. 

One  term  is  the  result  of  work  done  by  the  fluctuating  pressure  on 
the  vertical  motion  of  the  wall.  The  other  term  is  the  work  done  by 
the  viscous  forces.  These  terms  represent  energy  supplied  or 
absorbed  by  the  wall  and  must  be  compared  to  the  reduction  in  the 
energy  required  to  drive  the  main  flow  through  the  channel. 

In  Figure  16,  the  components  of  the  energy  balance  from  the  LES 
are  compared  to  the  results  of  the  linear  theory.  Three  cases  are 
represented  for  the  LES:  the  downstream- travel ing  wave  (case  1), 
and  two  upstream-traveling  waves  (cases  2  and  4).  The  LES  results 
are  generally  in  good  overall  agreement  with  the  linear  theory.  The 
pressure  work  term  is  particularly  well  predicted  by  the  theory, 
whereas  the  theory  estimates  the  viscous  work  to  be  higher  than  the 


LES  result. 


These  results  are  consistent  with  those  discussed 


previously  in  connection  with  Figure  7  (Section  3.4),  where  it  was 
shown  that  the  inclusion  of  the  nonlinearities  significantly  reduced 
the  gradient  of  the  turbulence  intensity  near  the  wall.  This  also 
results  in  a  reduction  in  the  viscous  work  term.  It  was  found  that 
the  intensity  gradient  was  further  reduced  in  the  LES. 

The  net  result  of  the  energy  balance  is  shown  in  Figure  17. 
Again,  the  results  of  the  linear  theory  are  reproduced  for 
comparison.  The  LES  results  seem  to  agree  with  the  linear  theory  in 
predicting  a  net  decrease  in  the  energy  requirement  for  both 
upstream-  and  downstream-travel ing  waves  in  the  low-speed  range. 
However,  the  calculation  of  the  net  energy  reduction  involves  a 
small  difference  between  quantities  which  are  of  similar  magnitudes, 
each  of  which  is  subject  to  some  uncertainty  due  to  numerical  error 
and  statistical  noise.  Although,  it  cannot  be  stated  unequivocally 
that  a  net  drag  reduction  is  produced,  the  fact  that  both  LES  and 
linear  theory  predict  such  a  reduction  makes  it  probable  that  a 
small  net  drag  reduction  is  actually  achieved. 


5.1.4  Pressure  Drag 


In  a  manner  similar 
linear  theory  in  Section 
calculated  in  the  LES  by 
the  fictitious  wall  over 


to  the  calculation  of  pressure  drag  by  the 
3.6  of  this  report,  the  pressure  drag  was 
integrating  the  pressure  fluctuations  at 
the  actual  wall,  i.e.. 


as 
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where  P  is  the  resolvable  portion  of  the  turbulent  wall-pressure 
fluctuation  from  the  LES  calculation.  The  results  of  this 
calculation  are  included  in  Table  1  as  compared  to  the  drag  of  the 
unperturbed  channel.  For  cases  1  and  4,  the  pressure  drag  predicted 
by  the  LES  is  somewhat  smaller  than  predicted  by  the  linear  theory. 


These  comparisons  are  only  a  qualitative  estimate  of  the 
channel  drag  since  there  is  no  actual  pressure  drag  with  the 
linearized  boundary  conditions.  The  interaction  of  the  pressure 
fluctuations  with  the  real  wall  cannot  be  evaluated  except  by  use  of 
a  calculation  scheme  that  employs  the  complete  nonlinear  boundary 
conditions.  This  needs  to  be  done  if  valid  results  are  to  be 
obtained . 


V-'  --'J 

4 


5 . 2  LES  Results  for  Compliant  Wall  Motion 

5.2.1  Effects  of  Turbulence  on  Wall  Motion 

Some  calculations  were  performed  for  compliant  walls  using  the 
coarse  mesh  LES  program  coupled  to  the  compliant  wall  Equations  (48) 
and  (49).  Several  cases  were  calculated  using  various  values  of  the 
elastic  constants,  kj,  C^,  and  and  membrane  tension  Nx  and  Nz 

(table  4).  The  first  case  used  values  given  by  Yang  and  Heller  as 


I 
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characteristic  of  a  Mylar  membrane  and  a  PVC  plastisol  substrate. 
The  non-dimensional  tension  forces  applied  to  the  membrane  were 
varied  from  0  to  10^  (corresponding  to  0  to  300  n/m).  The  material 
properties  from  Yang  and  Heller  (reference  19)  correspond  to  the 
following  non-dimensional  parameters 

1<3  =  k3T/(pu^)  =  1514 
C3  =  C3/(put)  =  .098 

D1  =  D1/(pu^3]  =  1.58 

Nx=  Nz  =  Nx/(pu^  6)  =  9082 


Figure  18  shows  the  time  history  of  the  rms  wall  deflection  (in  wall 
units)  and  the  rms  wall  velocity,  nt/uT  for  the  first  case  where 
Nx=  Nz  =  9082.  The  calculations  appear  to  have  reached  a 
statistically  steady  state.  However,  the  magnitude  of  the  wall 
deflection  is  very  small,  on  the  order  of  0.02  wall  units. 

As  the  prescribed  wall  motion  cases  had  a  nondimensional  surface 
deflection  for  the  prescribed  wall  shape  of  five  wall  units  and 
produced  a  small  effect,  this  deflection  will  produce  entirely 
negligible  effects.  As  should  be  expected  with  such  a  small  surface 
response,  no  effect  on  the  mean  pressure  gradient  was  observed. 
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Insight  into  the  nature  of  the  interaction  is  obtained  by 
examining  the  dispersion  relation  for  the  surface.  For  the  general 
viscoelastic  composite  material  described  by  Eg.  (47),  with  damping 


coefficient,  C3  =  0,  the  natural  frequency  is, 


0)  = 
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where  nx,  and  nz  =  0,1,2  .  .  .,  M  ,  =  Mj/  p~E  and  a  and  b  are  the 

length  and  width  of  the  membrane. 

For  a  simple  membrane  with  k3  =  0  and  negligible  stiffness,  the 
natural  frequency  becomes 


u 


“  -  t#'V^>2+  v^>2])1/2 


M. 


(54) 


In  the  present  case 


a  =  2b  =  2tt 
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The  characteristic  time  for  the  various  modes  can  be  found  from 


Tc  -  2 it/ oi 


For  N  /M  =  N  /M  =  2000 

X  1  Z  1 


T  =  2 TT  [ - -  -002 - =-]1/2 

V  (2"z> 


Values  of  Tc  for  the  first  6  modes  are  listed  in  Table  5. 


The  parameters  used  for  the  calculations  shown  in  Figure  18 
correspond  to  a  natural  frequency  of  a  spring-mass  system  which  is 
approximately 


-  TT 

u,  =  (-2)1/2-  45 

Ux 
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where  u*  is  the  dimensional  frequency  in  Hertz. 

An  estimate  can  be  made  as  follows  for  the  frequency  of  the 
turbulence  pressure  that  contains  the  maximum  energy.  From 
experimental  measurements  (reference  17)  the  convection  velocity  of 
the  turbulent  pressure  fluctuations  at  the  wall  is  estimated  to  be 
.  8Um«  A  nominal  value  for  the  mean  streamwise  extent  of  the  eddies 
producing  the  pressure  fluctuations  is  two  hundred  wall  units;  the 
corresponding  value  for  the  non-dimensional  frequency  is 

u  ~  300 

Anticipating  that  resonance  between  the  turbulent  pressure 
fluctuation  spectrum  and  the  vibration  characteristics  of  the  wall 
is  necessary  to  produce  significant  effects,  the  parameters  of  the 
surface  equations  were  adjusted  to  make  the  natural  frequency  of  the 
surface  more  closely  correspond  to  the  characteristic  frequency  of 
the  fluid  flow.  The  membrane  tensions  were  set  to  zero 
<NX  *  Nz  =  0),  the  viscoelastic  damping,  C3,  was  zero,  the  bending 
stiffness  Dj  was  increased  and  the  spring  constant,  k3,  was 
increased  to  raise  the  vibration  frequency  (case  2  in  Table  4).  The 
resulting  mean  rms  wall  deflection  was  n+  ~  .013,  indicating  no 
resonance  between  the  wall  and  the  flow,  and  as  before,  there  was  no 
noticeable  effect  on  the  mean  channel  pressure  gradient.  This 
indicates  that  the  estimate  of  the  required  resonant  frequency  was 
incorrect.  For  these  small  wall  deflections  the  linearized 


treatment  of  wall  boundary  conditions  is  accurate. 

Other  calculations  were  made  with  the  surface  equations 
containing  only  the  membrane  tension  terms.  That  is. 


k3  =  c3  =  D1  =  Q 


For  the  first  such  case  (case  3  in  table  4) 


N  N 

=  200 

M1  M1 

This  produced  a  wall  vibration  response  to  the  turbulent  pressure 
fluctuations  which  increased  continuously  (Figure  19),  showing  no 
tendency  to  stabilize.  The  calculation  was  halted  when  the 
amplitude  of  the  wall  deflections  reached  five  wall  units  due  to 
anticipated  errors  from  linearization  of  the  boundary  conditions. 
Increasing  N^/M^  an<^  ^z^l  t0  ^000  produced  a  wall  vibration 
which  increased  stepwise  (Figure  20).  Starting  from  no  vibration, 
the  wall  amplitude  increased  and  appeared  to  stabilize  at  an  rms 
value  of  two  wall  units.  After  about  600  calculation  steps  (a  time 
of  about  0.2256/ut),  the  amplitude  again  began  to  increase,  and  then 
stabilized  at  about  four  wall  units.  After  another  short  period  of 
stability,  the  amplitude  again  began  to  increase. 
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In  order  to  obtain  more  insight  into  the  nature  of  the 
interaction  the  pressure  fluctuations  and  the  wall  shapes  for  the 
various  cases  calculated  are  examined.  In  Figure  21,  the 
distribution  of  pressure  fluctuations  on  the  compliant  wall  is 
plotted  for  the  case  of  a  membrane  with  N^/M^  =  T^/M^  =  2000  at 
integration  step  600  (T  =  .225);  this  pressure  distribution  is 
typical.  It  is  quite  three-dimensional  and  contains  a  variety  of 
wavelengths.  By  contrast,  the  surface  shape  is  characterized  by 
longer  wavelength  displacements  and,  in  most  cases,  near  two- 
dimensionality  corresponding  to  the  natural  modes  of  vibration  of 
the  membrane.  Figure  22  shows  the  surface  displacement  for  the  Yang 
and  Heller  surface  which  has  an  rms  displacement  of  about 
.02  v/u^  (Figure  18).  Figure  23  shows  the  surface  shape  for  the 
modified  Yang  and  Heller  surface  (case  2,  table  4)  for  which  the 
damping,  C3,  was  zero,  the  membrane  tensions  were  zero  and  the 
spring  constant,  (<3,  was  adjusted  to  increase  the  natural  frequency 
of  the  elastic  material.  This  surface  shows  more  three- 
dimensionality  but  is  still  dominated  by  long  wavelengths  and  small 
amplitudes.  The  nearly  two-dimensional  surface  displacement  in 
Figure  22  corresponds  to  one  of  the  natural  modes  due  to  the 
dominant  membrane  terms  of  equation  (53),  whereas  the  more  three- 
dimensional  displacements  in  Figure  23  are  due  primarily  to  the 
spring  constant  term  l<3n  and  the  stiffness  term  DV^n.  Recall 
that  in  this  case,  the  membrane  stretching  terms  were  absent  (Nx  = 


These  results  indicate  the  need  for  information  about  the 
space-time  spectra  of  the  pressure  fluctuations  in  order  to  "tune" 
the  surface  parameters  to  achieve  resonance  with  the  flow.  Such 
information  would  be  very  difficult  to  obtain  experimentally,  but  it 
could  be  obtained  using  LES.  Such  calculations  are  suggested  as 
possible  future  work. 

The  largest  surface  displacements  were  obtained  for  the 

membrane  surfaces.  Figure  24  shows  the  membrane  surface  for  the 

case  with  TJ  /R  =  31  /H,  =  2000  at  five  different  times.  The 

x  1  z  1 

membrane  deformation  is  nearly  two-dimensional.  The  stepwise 
behavior  of  the  rms  amplitude  of  the  membrane  vibration  discussed 
previously  is  found  to  correspond  to  a  progressive  development  of 
different  modes  of  the  vibration. 

The  dispersion  relation  Eq.  (55)  shows  the  reason  for  the  two- 
dimensionality  of  the  membrane  surface.  First,  the  computational 
region  is  rectangular,  so  that  the  lowest  mode  is  dominated  by  the 
lowest  mode  of  the  longest  side,  which  is  in  the  streamwise 
direction.  For  integration  times  shorter  than  the  characteristic 
time  of  the  second  longitudinal  mode,  (see  table  4)  both  streamwise 
and  spanwise  waves  develop.  Then,  as  the  time  for  the  second 
longitudinal  mode  is  reached,  at  approximately  T  =  0.142,  the 
surface  develops  the  wave  shown  in  Figure  24  (a).  As  the 
integration  proceeds  further,  the  subharmonic  wave  begins  to  develop 
and  becomes  complete  as  the  integration  time  exceeds  T  =  0.28 
Figures  24  (b),  (c)  and  (d).  Further  integration  increases  the 
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amplitude  of  the  subharmonic  wave  while  the  higher  frequency  waves 
in  both  the  streamwise  and  the  spanwise  directions  are  damped, 
yielding  the  smooth,  two-dimensional  wave  shown  in  Figure  24  (e). 

Comparison  of  the  compliant-membrane  calculations  with  the 
prescribed-wall  results  reveals  that  while  the  compliant-wall  shape 
was  roughly  two-dimensional  and  sinusoidal  like  the  prescribed  wall, 
with  comparable  amplitude,  the  dominant  wavelength  produced  (the 
region  size)  was  much  longer  than  the  prescribed  value  (.05  x  region 
size).  The  prescribed  wall  motion  produced  a  discernible  drag 
reduction  because  it  had  the  proper  length  scales  to  interact  with 
the  energy-carrying  motions  of  the  flow.  The  compliant  walls  had 
wavelengths  too  long  to  interact  significantly  with  the 
turbulence.  To  match  both  frequency  and  wavelength  of  the 
turbulence  the  wall  wave  will  need  to  have  a  wavespeed  about  equal 
to  the  flow  speed.  This  is  the  essential  criterion  and  requires 
very  low  wave  speeds  for  solids.  So  the  wall  will  need  low  tension, 
high  density,  and  low  stiffness. 

5.2.2  Wall  Motion  Effects  on  Turbulence 

The  only  compliant  surfaces  for  which  noticeable  effects  were 

produced  in  the  flow  were  the  membranes.  Results  for  the  case 

with  N  /M  =  N  /M,  =  2000  are  shown  in  Figures  25-28.  In  Figure 

x  1  z  1  a 

25  the  turbulence  intensity  components,  (uj)',  (u^)',  (u^)'  ^or 
the  lower  and  upper  halves  of  the  channel  respectively  are  shown. 


As  for  the  prescribed  wall  case  shown  in  (Figure  9),  the  only 
significant  differences  occur  very  near  the  upper  wall,  where  the 
(up*  and  (up  '  increase  due  to  the  wall  motion.  Here,  the 
greatest  effect  is  seen  in  the  u^-component ,  while  the  u2-component 
has  only  a  slight  increase  in  intensity  due  to  the  wall  motion. 

In  Figure  26,  the  turbulent  shear  stress  is  shown.  In  contrast 
to  the  prescribed-wall  result,  there  is  no  noticeable  effect. 

Two-point  correlation  functions  as  defined  by  Eqs.  (50)  arl 
(51)  are  shown  in  Figure  27  at  the  same  distance  (x2  =  t  *990)  from 
each  boundary.  At  the  lower  boundary  [Figure  27(a)]  the 
longitudinal  correlations  are  typical  of  turbulent  flow  near  a  rigid 
boundary.  Near  the  upper  boundary  [Figure  27(b)],  the  longitudinal 
correlations  for  the  pressure  and  the  spanwise  velocity  component, 
U3,  are  unaffected  by  the  wall  motion,  while  the  u^  and  u2 
components  clearly  show  the  long  wavelength  of  the  boundary  shape. 

As  expected,  the  effect  did  not  propagate  very  far  into  the  flow. 

The  next  level  at  which  the  solution  was  examined  was  at  x^=  .941  ; 
no  effects  could  be  seen. 

In  the  spanwise  direction,  the  two-point  correlations  reveal 
the  two-dimensional  nature  of  the  wall  motion.  The  pressure  and  U3- 
velocity  correlations  are  unaffected  by  the  wall  motion,  but  the 
u^-  and  u2-components  are  highly  correlated  across  the  entire  span 
of  the  calculation. 

Energy  spectra  of  the  turbulence  are  shown  in  Figure  28  for  the 
same  two  locations  as  for  the  two-point  correlations. 


The  longitudinal  and  vertical  energy  components  show  a  significant 
increase  due  to  the  wall  motion  in  the  energy  contained  in  the  lower 
wavenumber  range. 


5.3  LES  Results  for  Other  Types  of  Wall  Conditions 

5.3.1  "Smart  Wall"  Experiments 

In  the  belief  that  the  optimum  drag  reduction  would  result  from 
a  wall  that  destroyed  all  of  the  turbulent  pressure  fluctuations, 
cases  which  imposed  such  a  boundary  condition  were  run.  In  the 
"standard"  problem,  the  three  velocity  components,  uj,  u2»  and 
U3,  are  specified  at  the  wall  and  the  pressure  is  computed.  In  the 
modified  calculation,  the  u^  and  U3  velocity  components  were 
specified  to  be  zero,  and  the  pressure  was  specified  to  be  constant 
at  the  wall  while  the  u2  velocity  component  was  computed.  This 
corresponds  to  a  porous  wall  with  a  plenum  behind  it.  Calculations 
resulted  in  a  higher  mean  pressure  gradient  indicating  that  the 
technique  actually  results  in  increased  drag. 

5.3.2  Distributed  Suction/Blowing 

A  boundary  condition  which  represents  a  distribution  of  suction 
and  blowing  through  a  wall  without  wall  motion  is 

V  u3  *  0 

U2  =  A  sin  (  ax  ^  tlt ) 


The  specific  formulation  used  here  has  a  simple  traveling  wave  form. 


Calculations  were  performed  for  this  formulation  both  with  LES 
and  the  mixing  length  analysis.  For  a  value  of  A  which  gives  the 
same  U2  amplitude  as  the  moving  walls,  the  Stokes  layer  Reynolds 
stress  is  smaller  than  when  the  accompanying  u^  boundary  condition 
is  present.  As  a  consequence,  the  effect  on  the  mean  pressure 
gradient  is  less  than  for  the  moving  wall.  Another  result  was  that 
the  change  in  the  mean  pressure  gradient  was  positive  for  waves 
moving  downstream  and  negative  for  waves  moving  upstream.  That  is, 
downstream-travel ing  waves  produced  a  drag  increase,  and  upstream¬ 
traveling  waves  produced  a  drag  reduction. 


5.3.3  Longitudinal  Grooves 


A  drag  reduction  technique  that  has  been  shown  to  be  effective 
experimentally  is  that  of  longitudinally-grooved  surfaces  having 
groove  depths  and  spacings  on  the  order  of  the  turbulent  wall-streak 
and  burst  dimensions  (reference  23). In  order  to  investigate  the 
applicability  of  LES,  two  cases  were  calculated  for  longitudinal 
grooves  of  sinusoidal  profile.  The  surface  shape  was  described  by 
Eq.  (11)  with  a  =  0  and  00  =  0 .  The  boundary  conditions  were  Eqs. 
(9)  and  (10).  The  values  used  for  the  other  parameters  were 


n  =5  and  6 
o 

0  =40  and  8 
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The  values  of  g  correspond  to  a  wavelength  (groove  spacing)  of  100 
and  500  wall  units,  respectively. 

The  results  were  that  the  short  wavelength  case  produced  a 
slight  increase  in  the  mean  channel  pressure  gradient,  while  the 
long  wavelength  case  showed  essentially  no  change  in  the  pressure 
gradient.  Experimental  results  (reference  23)  show  that  sharp-edged 
grooves  with  spacings  of  10-20  wall  units  are  required  to  produce  a 
significant  effect.  The  fact  that  no  effects  were  observed  in  the 
LES  with  smooth  grooves  is  consistent  with  those  results. 

6.0  LES  IN  GENERALIZED  COORDINATES 

It  was  shown  in  a  previous  section  that  the  use  of  linearized 
boundary  conditions  places  limitations  on  the  magnitude  of  the 
parameters  defining  the  wall  motion  and  especially,  on  the 
amplitude.  Furthermore,  the  values  required  to  produce  interesting 
effects  are  close  to  the  limits  of  validity  of  the  model.  In  this 
section,  a  formulation  of  the  equations  for  large-eddy  simulation  of 
the  channel  flow  using  the  exact  wall  boundary  conditions  is 
presented.  The  formulation  is  based  on  a  generalized  coordinate 
transformation  applied  to  the  Nav ier-Stokes  equations  and  then 
specialized  to  the  case  of  flow  in  a  channel  with  one  flexible 
wall.  Implementation  of  the  formulation  was  beyond  the  scope  of  the 
present  project,  but  is  required  if  LES  is  to  be  used  to  explore  the 
significant  part  of  the  parameter  space.  This  section  documents  the 


derivation  of  the  equations,  provides  some  analysis  of  potential 
sources  of  difficulty,  and  suggests  the  procedure  for  completing  the 
analysis . 


6 . 1  Governing  Equations 

The  complete  nonlinear  Navier-Stokes  equations,  with  velocities 

normalized  by  the  friction  velocity,  lengths  normalized  by  the 

2 

channel  half-width  and  pressure  normalized  by  puT  ,  can  be  written 


3D  =  _ 1 _ 3 _  f  3F  v 

3 1  3x  .  R  3x,  3x.  J 

1  e  k  k 


(  57) 


where 


6.2  Boundary  Conditions 


The  boundary  conditions  for  flow  in  a  channel  consist  of  the 
condition  of  no-slip  at  the  channel  walls  and  the  condition  that  the 
velocity  components  are  all  periodic  at  the  inflow  and  outflow 
boundaries  and  at  the  lateral  boundaries: 


ui  (xi ,t )  =  ®  f°r  =  1»2,3 
u^(x,h,x^»t)  =  u^ (x^  ,h,x^ ,t )  =  0 

3  h 

u  (x  ,h,x  ,t)  =  —  (x1,x3,t) 


u  .  ( x .  +  a,  X-  X-  t)  =  u  .  ( x ,  x  X  t) 

1  1  Z  ,  3  f  1  L  r  Z  ,  D  , 

u .  ( x ,  x„  x0  +  b,  t)  =  U.(x  x  X  t) 
1  1  t  Z  t  D  1  1  f  Z  ,  o  r 


where  h  is  the  instantaneous  distance  between  the  channel  walls  and 
a  and  b  are  the  length  and  width  of  the  computational  region  (note 
that  the  definition  of  x2  has  been  changed  so  that  x2  =  0  is  the 
lower  wall  and  x2  =  h  is  the  upper  wall). 


6.3  Coordinate  Transformation 


Define  a  general  transformation  of  the  form 


x  =  t 

=  f(Xj,t)  for  i,j  =  1  r  2 , 3 
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where  ^  is  defined  so  that  ^  =  0  is  the  channel  centerline  and  the 
channel  walls  are  at  ^  =  ±  !• 


The  chain  rule  yields 


3  3  3 
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2 

and 
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(63) 


The  inverse  transformation  is 


with  Equations  (62)  and  (63),  Equation  (57)  can  be  written 
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where  J  is  the  Jacobian  of  the  transformation 
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and  after  some  manipulation,  it  can  be  shown  that  equation  (66) 
becomes,  for  incompressible  flow. 
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6.4  Filterinc 


The  next  step  in  developing  the  flow  equations  for  large-eddy 
simulation  is  to  apply  a  filter  to  Equations  (69)  and  (70).  The 
filter  used  by  Moin  and  Kim  (reference  9)  will  be  applied  in  the 
transformed  coordinates.  This  is  equivalent  to  assuming  that  the 
small  scale  fluctuations  move  up  and  down  and  stretch  with  the 
moving  wall.  The  filter  takes  the  form 


f( VD  =  /J/v  f<VT>  i  n  g(c  -e‘)d^i 

k  =  l 


where,  for  i  =  1  and  3 
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and  for  1 
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6 . 5  Generalized  Equations 

It  is  assumed  that  the  wall  shape  consists  only  of  large  scale 
variations.  Applying  the  filter  to  Equations  (69)  and  (70)  yields 
the  dynamical  equations  of  the  large-scale  flow  field  in  the 
transformed  coordinates. 


The  continuity  equation  becomes: 
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(77! 


and  the  momentum  equations  are: 
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The  product  u.u.  contains  small-scale  components  which  must 

1.  K 

be  accounted  for.  Each  component  is  decomposed  into  a  filtered 
(large-scale)  part  and  an  unfiltered  (small-scale)  part,  eq.,  u^  = 
i  +  u'^.  Then  filtering  u^^  yields 


u .  u 

l 


k 


u  •  ui- 

l  k 


+  R-u 
lk 


(79 


where 
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The  subgrid  scale  stresses  are  defined  by 
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and  the  modified  pressure  is 


p  =  p  + 1  Rtt 


(82) 


The  filtered  momentum  equation  [(equation  (85)]  then  becomes 
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Numerical  solution  of  Equations  (77)  and  (83)  requires  careful 
consideration  of  the  conservation  properties;  these  are  discussed  in 
the  next  section. 

The  channel  flow  computer  program  developed  by  Moin  and  Kim  has 


been  thoroughly  tested  and  evaluated  and  has  been  found  to  be 
accurate,  and  stable  for  channels  with  flat  walls.  For  channels 


with  more  general  wall  conditions,  it  is  desirable  to  use  numerical 
methods  which  reduce  to  the  existing  method.  As  a  first  step  in 
developing  such  a  method,  define  transformed  velocity  components  as 


With  this 

(85) 

Using  the 


and  adding 
equation. 


definition,  equation  (72)  becomes 


30. 


(84) 


inverse  velocity  transformation, 


(86) 


appropriate  terms  to  both  sides  of  the  resulting 
equation  (83)  becomes,  after  some  manipulation 
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Equation  (87)  is  identical  in  form  to  equation  (5.1)  of 


(87) 


(88) 


reference  9.  The  terms  representing  the  velocity  components  and  the 


right-hand  side  of  the  equation  are  more  complicated  and  the 
turbulent  shear  stress  has  not  yet  been  completely  defined. 
The  boundary  conditions  are: 


0i=  0  at  ?2  =  -1;  i  =  1,2,3 


°1  =  °3  =  0  at  52  =  +1 

-  .  1  H2  1  3h  a.  .  _  .. 
°2  J  9x2  3t  C2  1 


Qi  =  periodic  in  $  and 


The  turbulent  shear  stress  in  equation  (88)  is  defined  in  the 
same  manner  as  in  reference  9.  That  is 


t..  =  2 v_( S .  .  -  <S.  >)  -  2v^  <S..> 
13  T  13  lj  T  13 


where  vT  and  v^.  are  the  small-scale  eddy  viscosities  which 
represent  the  action  of  the  unresolved  scales  of  motion  on  those 
that  are  resolved.  In  the  generalized  equations,  it  is  appropriate 
to  describe  the  shear  stress  in  terms  of  the  transformed  velocities 
and  coordinates.  Accordingly, 


6.6  Discretization  of  the  Generalized  Equations 


6.6.1  Transformation  for  a  Wavy  Wall 

It  has  been  pointed  out  by  several  previous  researchers  that 
when  studying  flows  over  wavy  walls  it  is  desirable  to  allow  the 
eddy  viscosity  and  turbulent  kinetic  energy  profiles  to  move  up  and 
down  with  the  wall.  Also,  since  measurements  are  most  easily  made 
with  a  probe  which  moves  up  and  down  with  the  wall  or  is  located  at 
a  fixed  point,  a  completely  curvilinear  coordinate  system  is  not 
practical.  A  coordinate  transformation  which  satisfies  these 
cons iderat ions  and  allows  the  calculations  to  be  performed  in  a 
channel  of  fixed  width  is 


Thus,  while  X2  varies  from  zero  to  hCx-^x^t),  ^  varies  from  -1 
to  +1 . 


6.6.2  Conservation  Properties 


It  is  generally  accepted  that,  in  the  numerical  treatment  of  a 
differential  conservation  law,  the  equations  must  be  in  strong 
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conservation  form  in  order  to  capture  discontinuities  properly  (cf. 
Lax  (reference  24),  Viviand  (reference  25),  and  Vinokur 
(reference  26)  and  others).  This  requirement  has  been  primarily 
enforced  for  compressible  flows  where  capturing  shock  waves  is 
important.  However,  these  properties  are  important  in 
incompressible  flows  as  well  (reference  27). 

In  the  incompressible  case,  because  of  the  uncoupling  of  the 
mass  and  momentum  conservation  equations  from  the  thermodynamic 
energy  equation,  the  momentum  equation  carries  a  second  important 
property;  conservation  of  kinetic  energy.  In  reference  9,  the 
momentum  equation  is  written  in  a  particular  form  which  is  found  to 
conserve  energy,  momentum,  and  circulation  with  virtually  any 
difference  scheme.  The  generalization  to  the  present  case,  is  not 
clear.  The  construction  of  an  energy-conserving  difference  scheme 
in  the  generalized  coordinate  problem  is  crucial,  and  should  be  the 
next  step  in  completing  the  analysis.  It  was  not  possible  due  to 
time  limitations  to  complete  that  step  in  this  effort. 

The  requirements  for  mass  conservation  and  the  time 
differencing  requirements  for  momentum  conservation  are  easily 
demonstrated.  Consider  equations  (85)  and  (87)  expressed  as 
difference  equations: 
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where  6/6C  is  a  finite  difference  approximation  to  3/9£  and  A 
represents  the  finite  difference  between  time  steps  n  and  n-rl , 
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Now,  conditions  which  will  ensure  conservation  in  a  uniform 
flow  field  are  examined.  If  the  flow  field  represented  by  the 
velocity  components,  uif  is  uniform,  then  Equation  (93)  with  (86) 
implies  the  following  difference  relation: 
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Similarly,  the  requirement  that  the  transformed  momentum  equation 
should  numerically  conserve  momentum  in  a  uniform  flow  leads  to  the 
following  difference  equation 
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(97) 


Now  the  term  H ^  contains  terms  identical  to  those  on  the  left-hand 
side  of  equation  (97).  Thus,  for  uniform  flows,  equation  (97) 
reduces  to 
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The  transformation  described  in  the  previous  section  yields  the 


following  difference  relations: 


Expanding  equation  (93)  gives 
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If  central  differencing  is  used,  then  the  transformation  yields 
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By  rearranging  the  right-hand-side  of  equation  (100)  and  using  the 
fact  that  <$h/65i  is  not  a  function  of  it  can  be  shown  that 


Thus,  if  central  differencing  is  used,  equation  (93)  will 
numerically  conserve  mass. 

The  condition  for  momentum  conservation  requires 
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For  the  transformation,  this  can  be  expressed  in  the  following 
difference  formulation: 
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Thus,  if  the  time  differencing  is  done  in  accordance  with 
equation  (104),  equation  (103)  will  be  satisfied  and  momentum  will 
be  conserved. 

6.6.3  Solution  Algorithm 

Since  equation  (87)  has  the  same  form  as  the  equations  used  in 
Reference  9,  the  numerical  method  used  in  that  reference  is 
suggested  as  the  one  to  use  first.  Thus,  the  pseudospectral  method 
would  be  used  in  the  5^  and  ^  directions,  and  central 
differences  would  be  used  in  the  direction.  Time  advancement 

would  be  made  using  a  semi-implicit  method  wherein  the  backward- 
difference  formula  would  be  used  for  the  time  derivative  of  the 
velocities,  the  Crank-N icholson  method  for  the  viscous  and  pressure 
terms  and  the  Adams-Bashf orth  method  for  the  remainder,  terms. 

7.  CONCLUDING  REMARKS 

A  program  for  large-eddy  simulation  of  turbulent  flow  in  a 
channel  was  adapted  to  deal  with  time-varying  wall  motion  through 
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the  use  of  linearized  boundary  conditions.  The  code  was  applied  to 
flow  in  channels  with  one  wall  having  prescribed  wavy  motion,  and  to 
channels  with  a  compliant  wall. 

A  mixing-length  analysis  employing  a  linearized  Stokes  layer 
theory  capable  of  qualitative  prediction  of  the  change  in  pressure 
gradient  was  developed  and  used  to  guide  the  selection  of  parameters 
for  LES  calculations.  The  results  indicate  that  a  small  drag- 
reduction  is  possible  under  certain  conditions  of  wall  motion.  The 
drag-reduction  is  the  result  of  the  creation  of  a  Stokes-like  layer 
whose  Reynolds  stress  is  opposite  to  that  of  the  undisturbed  flow. 
This  decreases  the  turbulence  production  near  the  wall  and  reduces 
the  pressure  gradient  needed  to  drive  the  flow. 

Predictions  using  both  the  linear  theory  and  LES  indicate  that 
much  of  the  reduction  in  the  energy  required  to  drive  the  flow 
through  the  channel  is  offset  by  the  work  required  to  drive  the  wall 
against  pressure  and  viscous  stresses.  The  net  effect  appears  to 
yield  a  small  decrease  in  the  energy  required  to  drive  the  flow. 

Estimates  of  the  drag  reduction  obtained  from  the  mixing  length 
analysis  are  generally  higher  than  those  obtained  from  LES. 

However,  the  mixing  length  analysis  provides  a  simple  first  estimate 
of  the  net  drag  reduction  due  to  a  wavy  wall. 

To  resolve  the  rapid  variation  of  the  Reynolds  stress  near  the 
moving  wall  a  fine  mesh  is  required.  For  this  reason,  the  mesh  was 
refined  beyond  that  used  in  the  plane  channel;  the  changes  were 


small . 


The  most  significant  effects  of  wall  motion  on  the  flow  are 
found  within  a  layer  whose  thickness  is  comparable  to  the  viscous 
sublayer  thickness.  To  validate  the  results  obtained  using  the 
linearized  conditions  requires  that  the  exact  boundary  conditions  be 
applied  at  the  actual  wall  location.  This  in  turn  requires 
transformation  to  a  time-dependent  coordinate  system  which  renders 
the  calculation  much  more  complicated. 

Calculations  for  compliant  elastic  walls  coupled  to  the  flow 
were  unsuccessful  in  producing  significant  effects  on  the  channel 
drag  for  material  properties  representative  of  existing  foam  or 
plastisol  materials.  The  only  surfaces  which  yielded  a  notable 
response  were  thin  membranes  under  slight  tension.  The  membranes 
developed  low  frequency,  long  wavelength  surface  shapes  which  did 
not  interact  with  the  flow  to  create  drag  reducing  effects. 
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Effects  of  a  Wall  with  Progressive  Waves 
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Figure  10.-  Turbulence  intensities  in  a  channel 
prescribed  sinusoidal  wall  motion  (case  1) . 
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Figure  13.-  Comparison  of  near-wall  F.eynolds  stress  from 
LES  and  nonlinear  Stokes  layer  solution. 


Wavenumber, 

(a)  Yu  s  55  measured  from  the  rigid  wall  (x~= 

W  fc 

Figure  15.-  Streamwise  energy  spectra  for  a  channel 
with  prescribed  motion  on  one  wall  (Case  1) . 
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Figure  18.-  Time  history  of  rms  wall  deflection  and  rms 
wall  velocity  for  first  compliant  wall  calculation. 
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Figure  26.-  Turbulent  shear  stress  distribution  in  a 
channel  with  one  wall  compliant. 
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